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ABSTRACT 

We present a search for periodic gravitational waves from the neutron star in the 
supernova remnant Cassiopeia A. The search coherently analyzes data in a 12-day in- 



■^^ Australian National University, Canberra, 0200, Australia 

^^The University of Melbourne, Parkville VIC 3010, Australia 

•"'Cardiff University, Cardiff, CF24 3 A A, United Kingdom 

■^^University of Salerno, 84084 Fisciano (Salerno), Italy 

^^Embry-Riddle Aeronautical University, Prescott, AZ 86301 USA 

^^The University of Sheffield, Sheffield SIO 2TN, United Kingdom 

^''inter-Umversity Centre for Astronomy and Astrophysics, Pune - 411007, India 

Southern University and A&M College, Baton Rouge, LA 70813, USA 
^^University of Minnesota, Minneapolis, MN 55455, USA 
^^Cahfornia Institute of Technology, Pasadena, CA 91125, USA 
^^Northwestern University, Evanston, IL 60208, USA 
^^University of Rochester, Rochester, NY 14627, USA 
*®The University of Texas at Austin, Austin, TX 78712, USA 
^''Eotvos University, ELTE 27 Budapest, Hungary 

**Rutherford Appleton Laboratory, HSIC, Chilton, Didcot, Oxon 0X11 OQX United Kingdom 

^^University of Adelaide, Adelaide, SA 5005, Australia 

®°Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain 

^^University of Southampton, Southampton, S017 IBJ, United Kingdom 

^^National Astronomical Observatory of Japan, Tokyo 181-8588, Japan 

^^Institute of Applied Physics, Nizhny Novgorod, 603950, Russia 

^''University of Strathclyde, Glasgow, Gl IXQ, United Kingdom 

^^Hobart and Wilham Smith Colleges, Geneva, NY 14456, USA 

^^University of Sannio at Benevento, 1-82100 Benevento, Italy 

^''Louisiana Tech University, Ruston, LA 71272, USA 

Andrews University, Berrien Springs, MI 49104 USA 
^^McNeese State University, Lake Charles, LA 70609 USA 
^"Cahfornia State University PuUerton, Pullerton CA 92831 USA 
^^Trinity University, San Antonio, TX 78212, USA 
^^Rochester Institute of Technology, Rochester, NY 14623, USA 
^■'Southeastern Louisiana University, Hammond, LA 70402, USA 



- 6 - 



terval taken from the fifth science run of the Laser Interferometer Gravitational- Wave 
Observatory. It searches gravitational wave frequencies from 100 to 300 Hz, and covers 
a wide range of first and second frequency derivatives appropriate for the age of the 
remnant and for different spin-down mechanisms. No gravitational wave signal was 
detected. Within the range of search frequencies, we set 95% confidence upper limits of 
0.7-1.2 X 10~^^ on the intrinsic gravitational wave strain, 0.4-4 x 10""^ on the equatorial 
ellipticity of the neutron star, and 0.005-0.14 on the amplitude of r-mode oscillations 
of the neutron star. These direct upper limits beat indirect limits derived from energy 
conservation and enter the range of theoretical predictions involving crystalline exotic 
matter or runaway r-modes. This paper is also the first gravitational-wave search to 
present upper limits on the r-mode amplitude. 

Subject headings: gravitational waves - stars: neutron - supernovae: individual (Cas- 
siopeia A) 



Introduction 



Using data from the Laser Interferometer Gravitational- Wave Observatory (LIGO; Abbott 



et al. 2009e), the LIGO Scientific Collaboration (LSC) and Virgo Collaboration have published 



searches for periodic gravitational waves from three astrophysically distinct types of rapidly rotating 



neutron star. Searches have targeted non-accreting pulsars (Abbott et al. 2004b, 2005b, 2007c 



2010), most notably the Crab pulsar (Abbott et al. 2008b, 2009b, 2010), using data from LIGO's 



first five science runs (designated S1-S5). Two searches have targeted the accreting neutron star in 
the low- mass x-ray binary Scorpius X-1, using data from S2 ( [Abbott et al]|2007aD and S4 ( [Abbott 
et al.[[20b7b ). Other searches have been broadband all-sky surveys for as-yet undiscovered neutron 



stars, using data from S2 ( [Abbott et al.[[2005a[ |2007a| ), S4 ( [Abbott et al.|[2008a[ [2009a[), and S5 
( [Abbott et ar[[2009c|d[ ). 



In this paper, we present the first directed search for periodic gravitational waves from a 
known, isolated, non-pulsing neutron star. The search targets the central compact object (CCO) 



in the supernova remnant Cassiopeia A (Gas A). The remnant is estimated to be 3.4^'^g'^ kpc 



distant (Reed et al. 1995), and to have been born in the year 1681 it 19 (Fesen et al. 2006b). It is 
the second-youngest known supernova remnant in the Galaxy, and the youngest with a confirmed 
CCO (Reynolds et al. 2008, De Luca[[2008 ). The remnant and CCO have been extensively studied 



through electromagnetic observations, which we summarize in Section 1.1 



There is compelling evidence that the Gas A CCO is a neutron star ( 


Pavlov et al. 


2000 


Chakrabarty et al. 


2001 


Gotthelf &; Halpern 


2008 


Pavlov &; Luna 


2009 


Ho Sz Heinke 


2009). We 



argue, in Section [1.2[ that its age and youth make it an interesting target for a search for periodic 
gravitational waves. Its youth means that it has not been covered by all-sky surveys for periodic 
gravitational waves, which focus on spin-down timescales much longer than the age of the remnant 
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( Abbott et al.||2005"al|2007a 2008a, 2009a|c|d ). We describe the first search for gravitational waves 
from Cas A in Section [2| present the results of the search in Section [3j and discuss them in Section |4j 
(We shall often abbreviate "the Cas A central compact object" to "Cas A.") 

Since the rotation frequency, and hence the gravitational wave frequency, of Cas A is unknown, 
we search for periodic gravitational waves with frequencies between 100 and 300 Hz. At these 
frequencies, where the strain noise of the LIGO detectors is lowest, the search is designed to beat 



indirect upper limits on gravitational radiation based on energy conservation (see Wette et al 
and Sections 1.2 and 2.3). The search found no credible signal (see Section [3|. In the 



2008 



absence of a detection, we present 95% confidence upper limits on a gravitational wave signal from 
Cas A, assuming its frequency is within the searched band. Upper limits are given for the intrinsic 
gravitational wave strain ho, the equatorial ellipticity e, and the r-mode amplitude a. 

Within the searched frequency band, the upper limits presented in this paper beat the indirect 



upper limits, as expected. Cas A is now one of only a handful of neutron stars (see Abbott et al. 



2010) where the most sensitive upper limits on gravitational radiation have been obtained using 



gravitational wave detectors such as LIGO. This paper is also the first gravitational-wave search 
to present upper limits on the r-mode amplitude. The best upper limits on e (a few times 10^^) 
and Q (a few times 10~^) are within the range of some theoretical predictions (see Section [4]). 



1.1. Electromagnetic Observations 



The Cas A CCO was first discovered as an x-ray point source in first-light images taken by the 



Chandra X-ray Observatory (Tananbaum 1999). It was subsequently identified in other satellite 



data dating back to the year 1979; the x-ray flux appears to have been constant since then (Pavlov 



et al. 2000). Optical and near-infrared searches have not found the CCO, and have all but ruled 



out the presence of an accretion disk from fallback, and a binary companion (Ryan et al. 2001 



Kaplan et al. 2001: Fesen et al. 2006a; Wang et al. 2007). The absence of the latter is puzzling. 



since the light-echo spectrum of the supernova, of type lib, implies that the progenitor was stripped 



of hydrogen by a companion (Young et al. 2006 Krause et al. 2008). One possibility is that the 



progenitor was the product of binary companions merging during a common envelope phase (Krause 



et al. 2008) 



The CCO appears to be a neutron star with a low surface magnetic field (an anti-magnetar). 



Blackbody fits to the x-ray spectrum by Pavlov et al. (2000) and Chakrabarty et al. (2001 ) implied 



a temperature too high, and an emitting area too small, to be consistent with emission from the 
whole surface of a cooling neutron star, or from the inner region of an accretion disk around a black 



hole. Pavlov et al. (2000) proposed a model of a strongly magnetized neutron star (a magnetar) 



with hot polar caps, which would, however, lead to x-ray pulsations which have not been observed 
(see below). Light echoes from the explosion have been interpreted as signs of a flare, reminiscent 



of a soft gamma repeater, occuring in the year 1953 (Krause et al. 2005). Such a flare, if single 
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could also have been due to a one-time phase transition ( Gotthelf &: Halpern||2008 ) . More recently, 
however, this interpretation of the echoes, and thus the evidence for the flare, has been discounted 



(Kim et al. 2008 Dwek & Arendt 2008). Recently, Ho & Heinke (2009) combined previous x-ray 



spectra (Hwang et al. 2004 Pavlov Sz Luna 2009), carefully adjusted for instrumental effects, and 
fitted them to various light-element atmosphere models. Their best fit was for nearly isotropic 
emission from the entire surface of a neutron star with a carbon atmosphere, low magnetic field, 
and mass and radius within normal ranges. 



The spin period of the neutron star is unknown. McLaughlin et al. (2001) searched for radio 



pulses (including possible binary orbital periods as short as a few hours) and found none at periods 
as short as 1-10 ms, depending on dispersion measure, making the CCO much more radio quiet 
than any known radio pulsar under 10^ years old. Searches for X-ray pulsations at periods as short 
as 2 ms ( |Chakrabarty et al.||200T| [Murray eral][2002{ [Mereghetti et"aL][2002| [Pavlov k Luna||2009 ) 
have produced at best a marginal candidate for a period at 12 ms, which has not been confirmed. 



No pulsar wind nebula has been detected (Hwang et al. 2004, Pavlov &; Luna 2009). 



1.2. Motivation for a Gravitational Wave Search 



If the CCO in Cas A is an anti-magnetar, it may be spinning fast enough to emit periodic 
gravitational waves above 100 Hz, where LIGO is most sensitive. 

To date, seven supernova remnant CCOs are known: three have observed spin periods, and 



a fourth has an observed periodicity that may be due to binary orbital motion (De Luca 2008 



Gotthelf & Halpern 2009). The fastest-spinning CCOs have rotation periods of ~ 100 ms, which 
are much slower than the longest rotation periods covered by this search (20 ms for gravitational 
waves from a non-axisymmetric distortion, and 13 ms for r- modes; see Section [s]). Only one of 
the CCOs has a measurable spin-down ( Halpern Gotthelf[[2010 ) and the others have tight upper 
limits. This indicates that the CCOs have magnetic fields much less than those of typical radio 
pulsars, and rotation periods that have not changed significantly since birth. If Cas A spins as 
slowly and constantly as these objects, it is not detectable by LIGO. It is difficult, however, to 
definitively extrapolate the general properties of CCOs from a sample of three or four. 

Young neutron stars, such as Cas A, may be the most likely to retain non-axisymmetries from 
the violent circumstances of their births. For example, the formation of the crust during an epoch 
of perturbations such as r-modes ( Lindblom et al.[[2000 Wu et alT][2001 ) could lead to an irregular 
shape. Gravitational radiation may also be generated by the continuing non-axisymmetry of the 



r-modes themselves (Owen et al. 1998), which may last for up to thousands of years (Arras et al. 



2003) depending on the composition and viscosity of the star. 

Indirect upper limits on gravitational wave emission from Cas A can be estimated using a 



(Wette et al. 


2008 


Owen 


2010 



star was born spinning at least ~ 20% more rapidly than it is now, and the spin-down evolution 
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has been dominated by the emission of gravitational waves, the unknown spin frequency and fre- 
quency derivative can be ehminated in favor of the known age to place a rough upper limit on the 
gravitational wave emission. For Cas A, the indirect limit on the intrinsic gravitational wave strain 



ho (see Section 2.2) is 



ho < 1.2 X 10-24n-ir. 



1/2.1/2 

/I 



3.4 '300 -"45 



(1) 



where 1)3.4 is the distance to Cas A in units of 3.4 kpc, T300 is its age in units of 300 yr, and is 
its princip 
frequency. 



its principal moment of inertia in units of 10^^ g cm^. The indirect limit on ho is independent of 



The choice of a fiducial age of 300 yr for Cas A, at the young end of the range estimated by 



Fesen et al. (2006b), is conservative in that it gives a larger search parameter space (see Section 2.3). 



It also raises the indirect limit by ~ 10%, a small effect compared to the uncertainties in the distance 
(of order 10%), and the principal moment of inertia, which may be up to 3 times higher than its 



fiducial value (see Abbott et al. 2007c). The uncertainties in the direct upper limits presented in 



Section [3] are on the order of 10-15%, due to uncertainties in the calibration of the LIGO detectors 
(of order 10%; see Section [2.1[ ), and systematic uncertainties in the search pipeline (of order 5%; 
see Section 2.2). Equation ([T]) assumes gravitational waves from a mass quadrupole; the equivalent 



limit for r- modes is higher by tens of percent (Owen 2010) 



The indirect limit on ho may be converted into an indirect limit on the equatorial ellipticity 



f < 3 q X 10-4^-1/2.-1/2 .-2 



(2) 



where /loo is the gravitational wave frequency in units of 100 Hz. Due to the uncertainties in 
neutron star parameters mentioned above, this limit on e is overall uncertain by roughly a factor 2. 



The indirect limit on ho also implies an indirect limit on gravitational waves from r-modes (Owen 



2010), in terms of their amplitude 



a 



< O.Mrg-yV: 



3 

100 ' 



(3) 



where the uncertainty due to the moment of inertia and other properties of the star, while more 
complicated, is roughly a factor of 2-3. If a varies over the observation time, this limit applies to 
the rms value of a over time. See |Owen| ( |2010D for precise definitions of a and e, translations to 
other quantities used in the literature (including ho), and more discussion of uncertainties. 



Wette et al. (2008) showed that a search for Cas A of 12 days of LIGO S5 data in the band 



100-300 Hz is feasible and can expect to beat the indirect limits. The Cas A indirect limits on 
ho and e are comparable to the spin-down limits for the Crab pulsar, which have been beaten 



by searches of LIGO S5 data (Abbott et al. 2008b, 2010). The indirect limits on a are lower and 



therefore more interesting than for the Crab if Cas A is emitting at 100-300 Hz. The indirect limits 
are also comparable to the best upper limits achieved by all-sky searches for periodic gravitational 
waves (Abbott et al. 2009c|d ). The all-sky searches covered longer spin-down timescales and thus 
were not sensitive to Cas A if it is emitting gravitational waves near the indirect limit. 
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2. Search Pipeline 

This section describes the gravitational wave search for Cas A. The search pipeline consists 
of: the selection of data for the search (Section 2.1), the analysis method (Section |2.2[ ), the search 
parameter space (Section 2.3), the template bank used to compute the search (Section 2.4), post- 
processing of the resuhs (Section 2.5), examination of the significance of the largest value of the 
detection statistic returned by the search (Section 2.6 ), and, in the event of no detection, calculation 
of the upper limits (Section 2.7). 



2.1. Data Selection from the LIGO S5 Run 

LIGO is a network of three interferometric detectors: a 4-km arm-length detector in Livingston, 
Louisiana (LI) and two detectors of 4-km (HI) and 2-km (H2) arm lengths co-located in Hanford, 
Washington. 



The S5 science run (Abbott et al. 2009e) is LIGO's fifth and most recently completed science 
run. It commenced at 2005 November 4, 16:00 UTC at Hanford, and at 2005 November 14, 
16:00 UTC at Livingston; it ended at 2007 October 1, 00:00 UTC. The S5 run collected over one 
year of science data coincident among all three detectors, with an overall triple-coincidence duty 
cycle of 54%. Interruptions caused by environmental disturbances, as well as scheduled breaks 
for maintenance and commissioning of equipment, accounted for the downtime. During S5, the 
detectors were operating at very near their design sensitivities. The strain noise of the two 4-km 
detectors was on average less than 3 x 10"'^'^ Hz~^/^ at their most sensitive frequencies (around 
140 Hz) and less than 5 x 10"^^ Hz~^/^ over 100-300 Hz, and generally improved (as did the duty 
cycle) over the course of the run. 

This search uses science data from only the LI and HI detectors. Data from the H2 detector 
is less sensitive, but carries the same computational cost to search. A small percentage of the 
acquired science data is excluded by data quality controls, which identify times when the data is 
known to be unsuitable for analysis. This includes, e.g., data taken when the output photodiodes 
of an interferometer were saturated, when the calibration of the data was ill-defined, when high 
winds were measured at the Hanford observatory, and 30 seconds before an interferometer lost lock. 
The remaining science data is calibrated ( Abbott et al.||2010 2009e; Abadie et al.||2010 ) to produce 



a (discontinuous) time series of gravitational wave strain, h{t). The time series is then broken 
into 30-minute segments. Because not every continuous section of h{t) is an integer multiple of 30 
minutes in length, some science data is discarded in this process. Finally, each 30-minute segment is 
high-pass filtered above 40 Hz and Fourier transformecQ to form Short Fourier Transforms (SFTs) 
of h{t). The SFTs are the input data to the search pipeline. 



^ Prior to the Fourier transform, each segment is multiphed by a nearly-square Tukey window to mitigate transients 
at the start and end of the segment. The loss in signal power due to the windowing is on the order of 0.1%. 
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The maximum uncertainties in the caHbration of h{t) are 10.4% in ampHtude and 4°. 5 in phase 
for HI, and 14.4% in amphtude and 4°. 2 in phase for LI; all uncertainties are constant in time 



to within 1% (Abadie et al. 2010). The uncertainties in the calibration amplitude contribute to 



the overall systematic uncertainty in the upper limits presented in Section [3] The analysis method 



used in this search (see Section 2.2) is sensitive only to the relative phase uncertainty between 
detectors. Even for a worst-case relative phase uncertainty of ~ 10°, the resulting difference in the 
signal frequency between detectors would still be much less than the mismatch allowed for by the 



search template bank (see Section 2.4). Thus, the phase uncertainties do not affect the results of 
this search. 



The search is restricted to a data set spanning a maximum of 12 days (Wette et al. 2008). 
The computational cost of a coherent search for Cas A scales with the 7th power of the timespan 
between the first and last timestamp of the analysed SFTs, while the sensitivity scales only with 
the square root of the observation time. Increasing the timespan from 12 days would therefore 
rapidly increase the computational cost, for a negligible improvement in sensitivity. To select the 
12-day data set, we compute a figure of merit for each possible data set spanning 12 days, chosen 
from the available S5 SFTs. The figure of merit, which is proportional to the estimated power 
signal-to-noise ratio, is given by "^j^ f[Sh{f)]~^ , where Sh is the strain noise power spectral density 
at frequency / in the SFT numbered k, and the summation is over all SFTs in the data set and 
over the frequency band 100-300 Hz. We select the data set with the maximum value of the figure 
of merit. 

At the time this search was conducted, SFTs were available from the beginning of S5 until 
2007 April 19 UTC. In this period, ~ 9% of the HI science data and ~ 13% of the LI science 
data is excluded, either due to data quality vetoes or due to the segmentation of h{t) during SFT 
generation. The 12-day data set selected for the search begins at 2007 March 20, 20:56:37 UTC 
and ends at 2007 April 1, 20:50:04 UTC. It contains a total of 934 SFTs (445 from LI and 489 
from HI) and an average of 9.7 days of data from each detector. 



2.2. Analysis Method 

The data is searched using the J-"-statistic, a coherent matched filtering technique used to search 



for periodic gravitational waves using multiple detectors (Jaranowski et al. 1998 Cutler & Schutz 



2005). Matched filtering requires an accurate model, or template, of the signal. The template 
models the response of a detector to the two polarizations ("-I-" and "x") of the gravitational 
wave signal emitted by a rotating neutron star. In addition to the sky position, the barycentered 
gravitational wave frequency, and its derivatives, the signal template has four parameters related 
to amplitude and polarization: the intrinsic strain ho, initial phase constant (pQ, inclination angle 
L of the star's rotation axis to the line of sight, and polarization angle ip. The J-"-statistic is the 
logarithm of the likelihood ratio analytically maximized over these unknown parameters. The value 
of the J^-statistic is usually quoted as 2J^. 
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In transforming a signal from the reference frame of the star to that of the detector, the 
intrinsic frequency of the source is modulated by Doppler effects due to the sidereal and orbital 
motion of Earth with respect to the source. The transformation requires the right ascension a and 
declination 6 of the star, which for Cas A are known to high precision: a = 

and 6 = 58°48'42".51 ± 0".4 (iFesen et al.||2006ap. The sky resolution of the J"-statistic is, for a 



data set spanning 12 days, much coarser than the measured uncertainties in a and S (Whitbeck 



2006), and so no search over sky position is required. The search is conducted over the remaining 
unknown function in the signal template: the instantaneous frequency of the source as observed at 
the solar system barycenter, f{t). This is modeled as 



f{t)^f + f{t 



to) + lf{t 



to)' 



(4) 



where the initial frequency /, first spin-down /, and second spin-down / (all evaluated at the start 
time of the data set, t = to) constitute the search parameters. In previous searches for periodic 
gravitational waves, it has not been necessary to include a second spin-down, but one is required 



for this search due to the young age of Cas A (Wette et al. 2008). 



The signal template does not allow for the possibility that Cas A glitched during the 12 days 
spanned by the data. On the other hand, even the most frequently glitching pulsar does so only a 
few times per year, and glitch population statistics do not clearly indicate that the youth of Cas A 



is sure to mean more frequent glitches (Yuan et al. 2010). The worst case would be a glitch at the 



midpoint of the 12 days span, since in other cases the template would pick up the longer of the 
pre- and post-glitch coherent stretches. Amplitude signal-to-noise accumulates as the square root 
of observation time, so the worst loss would be a factor of ~ 2. 

The J^-statistic implicitly assumes a uniform prior on the polarization angle ip, and therefore 
the standard data analysis and upper limit procedures remain valid for gravitational waves signals 
from r-modes, even though the signal template used in this search assumes a mass quadrupole 



(Owen 2010). 



The Cas A search uses the ComputeFStatistic_v2 implementation of the J-"-statistic, which 
is available as part of the LALSuite software packagej^ Values of the J^-statistic returned by 
ConiputeFStatistic_v2 have an uncertainty of up to 5%, due to practical computation issues and 



optimizations (see Prix 2010). 



2.3. Parameter Space 

The range of the gravitational wave frequency, 100 Hz < / < 300 Hz, is chosen based on the 



estimate in Wette et al. (2008) of the frequency band over which a search of LIGO S5 data could 



■^The version of the software used in the search is tagged with the identifier S5CasASearch. See https : //www . 
Isc-group .phys .uwm. edu/daswg/projects/lalsuite .html 




Fig. 1. — Visualization of the Cas A search parameter space. The black-outhned shapes are sHces 
of the /-/ parameter space at the following fixed values of /: (back to front) 100, 150, 200, 250, 
and 300 Hz. 
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beat the indirect limit on /iq at reasonable computational cost. 



The ranges of the spin-down parameters / and / are chosen to be (Wette et al. 2008) 



/ 



(minn — 1)t 



/ 



(maxn — 1)t 



and 



(minn)/^ </ < ^^^^^^■f'^ 



(5) 



(6) 



/ - f 

respectively, where the braking index n is defined below. The age of Cas A, r, is chosen to be 



300 yr (as discussed in Section 1.2). Note that the range of / depends on /, and the range of / 



depends on both / and /. The resulting shape of the three-dimensional parameter space of /, /, 
and / is depicted in Figure [l] 

The choices of / and / ranges are motivated by the desire to cover a wide range of astro- 
physically motivated possibilities, expressed in terms of the braking index. In equation ([g]), we use 
the definition of the instantaneous braking index n = fff~'^; in equation ([s]), the angled brackets 
denote the average value of n over the lifetime of Cas A. If the dominant emission mechanism 
driving the spin-down of Cas A has changed over its lifetime, these two quantities will be differ- 
ent. Therefore, in the spirit of trying to cover the broadest imaginable parameter space, we do 
not constrain the instantaneous and averaged values of n to be the same. Instead, we search over 
values of n (both instantaneous and averaged) between 2 and 7. This range covers the possibilities 
that Cas A is spinning down primarily due to magnetic dipole radiation (n = 3) or to gravitational 
waves generated by a mass quadrupole (n = 5) or by constant-a r- modes (n = 7). This range 
also covers the braking indices of nearly all known pulsars, which are typically between 2 and 3 



(Livingstone et al. 2006). The exception is the Vela pulsar, for which n ~ 1.4 (Lyne et al. 1996) 



Extending n to lower values would dramatically increase the computational cost of the search. 



2.4. Computation 



The search consists of computing the J^-statistic over a finite bank of templates, whose pa- 
rameters are given by points within the search parameter space. The search uses a template bank 



generation algorithm (Wette 2009) which locates the parameter space points at the vertices of a 



body-centered cubic lattice. This minimizes the number of points per unit volume required to cover 
the parameter space (see, e.g., Conway &: Sloane|1988 Prix|2007b ). The points are spaced using the 
J-'-statistic parameter space metric ( |Whitbeck|[2006l |Prix|[20b7a| ) to ensure that the maximum ex- 
pected fractional loss in signal-to-noise ratio (known as the mismatch) will never exceed 20%. Due 
to strong correlations between the frequency and spin-down parameters, and the irregular shape of 
the parameter space, it was necessary to place additional templates outside of the parameter space 
to fully cover the parameter space boundaries, particularly for /. As a result, the number of tem- 
plates searched, ~ 7 x 10^^, is an order of magnitude larger than that estimated in Wette et al. 
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(2008). This resulted in a greater computational cost but does not greatly affect the sensitivity of 



the search, which is only weakly dependent on the number of templates (see Section 2.6). 



The search is divided into ~ 21500 independent computational jobs by partitioning the range 
of / into small bands, each with approximately equal numbers of templates, resulting in band 
widths of 1.8 mHz (at / = 100 Hz) to 20 mHz (at / = 300 Hz). The search completed in ~ 3.5 
days on ~ 5000 cores of the ATLAS computer cluster at the Max Planck Institute for Gravitational 
Physics (Albert Einstein Institute) in Hanover, Germany. From each search job, only the 0.01% of 
templates with the largest values of 2T were recorded. 



2.5. Post-Processing 

Data from the LIGO detectors are known to contain stationary or nearly stationary spectral 
lines originating from instrumental and environmental noise. The J-"-statistic depends upon a robust 
estimator of the power spectral density of the noise. The implementation used in this search (see 



Section 2.2 ) uses a spectral running median with a window size of 50 SFT bins, or 27.8 mHz (Abbott 



Section 


2.2) 


et al. 


2007a 



et al. 2007a). Lines that are narrower than the median window size will remain in the data, and 



may result in spuriously large values of the J^-statistic. 

To identify frequency bands where the ^-"-statistic may have been contaminated in this manner, 
we search for prominent narrow lines in power spectra of the searched HI and LI data. This 
procedure identified eleven frequency bands, which are listed in Table [T} along with the largest 
value of the J-'-statistic found in each band. Templates whose instantaneous frequency /(t) at 
any time falls within any of these bands are excluded from the remainder of the search pipeline. 
Approximately 2 x 10® templates (a fraction ~ 3 x 10~® of the total number of templates) are 
excluded in this manner. Four of the bands contain a largest value of 2T in the range 51-59, which 



would not be regarded as statistically significant gravitational wave candidates (see Section 2.6). 



At certain times during S5, ten simulated periodic gravitational wave signals were injected 
into the LIGO detectors at the hardware level by mechanically oscillating the detector mirrors 



(Abbott et al. 2008a). Four of the injections had frequencies within the search band. The Cas A 
search data set was not selected with regard to times when the hardware injections were active. 
As a result, < 1 day of the searched data (~ 4.2% of the HI data, and ~ 9.2% of the LI data) 
contain the injections, and the effective intrinsic strains of the injections are reduced by a factor of 
~ 10. After accounting for this reduction, two of the injections have strains of /iq ^ 10~^^, and are 
undetectable by this search. The remaining two injections (designated nos. 3 and 8) have reduced 
strains of /iq ~ 1-63 x 10^^'^ and 1.59 x 10~^^ respectively. These injections are found by the Cas A 



search at frequencies consistent with the parameters of the injected signals (Wette 2009). Neither 
injection is at the sky position of Cas A, but are nevertheless detected due to their strength, and the 
poor sky localization of short-duration (i.e. less than 1 day) signals arising from global correlations 
in the signal parameter space dPrix Itoh | [2005l |Pletsch|[2008| ) 
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The remaining nine non-injection bands are ruled out as gravitational wave candidates for the 
following reasons. Two contain harmonics of the 60 Hz power mains frequency, and one contains a 
harmonic of the 16 Hz data acquisition buffering frequency (Abbott et al. 2004a). The remaining 
six bands contain narrow instrumental lines which occur only in the LI detector. We would expect 
a convincing gravitational wave candidate to be seen in both detectors. Four of these six lines 
(at ~ 139.2 Hz, 144.7 Hz, 185.6 Hz, and 193.0 Hz) are definitively identified with environmental 
noise, by correlating the gravitational wave data with data from environmental monitors (e.g. 
magnetometers, accelerometers and microphones). The sources of the remaining two lines, at 
~ 139.5 Hz and at ~ 209.2 Hz, are not conclusively identified. The value of the J^-statistic 
associated with the ~ 209.2 Hz line, 2J- = 54, is smaller than the largest value of 2J- expected from 
the search in the absence of a gravitational wave signal (see next section), and is therefore not a 
candidate for a gravitational wave signal from Cas A. The value of the ^-"-statistic associated with 
the ~ 139.5 Hz line, 2T = 72, is also not statistically significant; as may be deduced from Figure[2j 
there is a ~ 5% probability that the search would return a higher largest value of 2T without a 
gravitational wave signal being present in the data. 



2.6. Significance of the Largest 2T 

The largest value of the ^-"-statistic returned by the search after post-processing is denoted 2T* . 
This is our most promising candidate for a gravitational wave signal from Cas A. The probability 
density of 2J^*, under the assumption that no gravitational wave signal from Cas A is present in 
the searched data, is given by 



p(2^^) = 7Vp(xl;2^*) 



2T^ ^-1 

d(2^)p(xl;2^) 







(7) 



where A'" is the number of searched templates, and p(x4;2J^) denotes the probability density of 
a central distribution with 4 degrees of freedom (i.e. the distribution of 2T in the absence 
of any signal, and assuming Gaussian noise). If the value of 2J^* returned by the Cas A search 
(see Section [3]) is within the range of probable values expected from p(2J^*), it is not statistically 
significant. If, on the other hand, the value of 2T*' is extremely unlikely to be drawn from p{2J^*), 
the candidate signal is worthy of further investigation. 

Equation Q assumes that each of the values of 2J^ are statistically independent. This is not 
strictly true, however, as we expect the 2J^ values of neighbouring templates to be correlated due to 
the close template spacing. Therefore, in equation ([T]), N should be substituted with the number of 
statistically independent templates, Ni < N. An empirical estimate of the number of statistically 
independent templates found that Ni ~ 0.88 A for this search ( Wette|2009 ). A reduction in N shifts 



the distribution of p(2J^*) towards lower values of 2J^, and thus increases the statistical significance 
of a candidate signal. For a reduction of N to 0.88A, however, the shift is negligible, and can 
be ignored; indeed. Figure [2] plots equation ([T]) with Ni = N. If, taking an extreme example. 
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Ni = O.IA^, the position of the maximum of p[2T*) would be shifted from ~ 66 to ~ 62 (see 
Figure [2]) , but the significance of the found by this search (see Section [3| would not be greatly 
increased. Therefore, the uncertainty in the number of statistically independent templates does 
not alter the significance of 2J^*, nor the conclusions reached in Section [3} The determination the 
number of statistically independent templates from first principles is an interesting area for further 
investigation. 



2.7. Upper Limits 



If, after examining 2T* (see Section 2.6), we conclude that no gravitational signal has been 
detected, we proceed to set 95% confidence upper limits on the intrinsic strain /iq, the ellipticity 
e, and the r-mode amplitude a. The upper limits are determined using Monte Carlo injections, 



following the procedure described in Abbott et al. (2007a). The search frequency band is first 



partitioned into 400 sub-bands of width 0.5 Hz; an upper limit is set separately for each sub- 
band. The choice of 0.5 Hz is small enough that the time-averaged noise fioor of the detectors 
is approximately constant over each band, and is large enough to keep the computational cost 
reasonable. We denote the largest values of 2T found by the full Cas A search in each sub-band 
by 2T*_^] they serve as the false alarm thresholds in each band. For each sub-band, a population 
of 5000-6000 periodic gravitational wave signals, with fixed and all other parameters randomly 
chosen, are each in turn injected into the search data set and then searched for using the same 
analysis method used in the full Cas A search. We record the largest value of 2J^ found by the 
searches for each of the injected signals. The fraction of these 2T values that are greater than 2T^j^ 
gives the confidence corresponding to the fixed value of /iq. 

To expedite the injection procedure, we use an analytic model of the distribution of the popu- 
lation of injected signals to provide an educated guess at the value of /iq required for 95% confidence 
( Wette|2009 ). We then perform the injection procedure once, at that /iq; to check that the recovered 



fraction really is 95%. We find that the analytic model slightly overestimates the required /iq, and 
so the 95% confidence upper limits presented in Section [3] are conservative. 



3. Results 



After post-processing, including the exclusion of the 11 frequency bands discussed in section 2.5 
the largest remaining value of the J^-statistic returned by the search is 2T* ~ 65. Figure [2] compares 
2T* to its expected theoretical distribution, under the assumption that the searched data contains 



no gravitational wave signal from Cas A (see Section 2.6). It is clear that 2T* is consistent with this 
distribution, and is therefore not statistically significant. We therefore conclude that the searched 
data does not contain any plausible gravitational wave signal from Cas A. This conclusion is not 
significantly influenced by the uncertainty in the number of statistically independent templates (see 
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Table 1: Frequency bands (column 1) identified during post-processing as containing spuriously 
large values of the ^-"-statistic (the maximum of which are given in column 2), and a brief description 
of their origin (column 3). See the text for details. 



Frequency Band 


2-7^ max 


Origin 


108.860 ±0.018 


90 


Pulsar hardware injection no. 3 


119.877 ±0.019 


72 


Sideband of 60 Hz harmonic 


128.000 ±0.017 


56 


16 Hz harmonic 


139.225 ± 0.058 


70 


Ll-only line 


139.510 ±0.017 


72 


Ll-only line 


144.751 ± 0.056 


110 


Ll-only line 


179.812 ±0.018 


51 


Sideband of 60 Hz harmonic 


185.630 ±0.055 


59 


Ll-only line 


193.005 ± 0.046 


66 


Ll-only line 


193.391 ±0.018 


73 


Pulsar hardware injection no. 8 


209.265 ±0.017 


54 


Ll-only line 




55 60 65 70 75 

Fig. 2. — The solid line is the probability density of the largest value of the J-"- statistic, 2J^*, under 
the assumption that no gravitational wave signal from Cas A is present in the searched data [see 
equation ([T])]. The dashed line denotes the value of returned by the Cas A search. 
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Section 2.6). 

Upper limits at 95% statistical confidence on ho, e, and a are plotted, alongside their respective 
indirect limits, in Figures[3j|4j and[5]respectively. Systematic uncertainties in the direct upper limits 
on ho are of order 10-15% (see Sections |2. 1| and |2.2[ ) . Systematic uncertainties in all three indirect 
limits, and in the direct upper limits on e and a inferred from ho, are roughly a factor 2-3 for 



a and 2 for the others (see Section 1.2). As expected, the upper limits beat the indirect limits 
over the gravitational wave frequencies 100-300 Hz. The equivalent ranges of rotation periods are 
6.7-20 ms for the upper limits on ho and e (which assume the gravitational wave frequency to be 
twice the rotation frequency), and 4.4-13 ms for the upper limits on a (which assume a 4/3 ratio 
of gravitational wave frequency to rotation frequency). Upper limits within 2 Hz of harmonics of 
the 60 Hz power mains frequency are severely degraded due to noise, and are excluded from the 
figures. Between 230-240 Hz, and around 280 Hz, the upper limits are degraded by disturbances 
in the broadband noise of the HI detector. 



Discussion 



The tightest upper limit on ho (Figure [3j) is ~ 7 x 10 at ~ 150 Hz, in the region where 
the LIGO detectors are at their most sensitive. The search improved slightly upon the expected 



upper limits on ho estimated by Wette et al. (2008), due to slightly better detector sensitivity and 



duty cycle in the selected data set. Therefore, the search could have beaten the indirect limits 
at frequencies slightly outside of the 100-300 Hz band, although the computational cost increases 
rapidly at higher frequencies, and the noise floor of the detectors rises steeply at lower frequencies. 

The upper limits on e (Figure |4]) range from ~ 4 x 10~^ at 100 Hz (20 ms rotation period) 
to ~ 4 X 10"^ at 300 Hz (6.7 ms), assuming the canonical parameters of Eq. Q. The upper 
limits are higher than the maximum e of a few times 10~^ predicted for normal neutron stars. 



even with recent results indicating a high breaking strain of the crust ( Horowitz &; Kadau 2009 ) . 
Ellipticities of a few times 10~^ are within the range of predictions ( Owen 



2005 



Li^|2007| 



Haskell 



et al. 2007 Knippel & Sedrakian 2009) for various forms of crystalline quark matter (Xu 2003 



Mannarelli et al.||2007[ ). Robust hybrid models ( Glendenning||1992 ) could sustain ellipticities up to 
about 1 X 10~^ (scaled from Owen||2005 ) if the breaking strain of Horowitz & Kadau (2009) is valid 
for the mixed phase of matter. Ellipticities comparable to our upper limits could also be sustained 
by internal magnetic fields of order 10^^ G, depending on the field configuration, equation of state. 



and superconductivity of the star (Cutler 2002 Haskell et al. 2008; Akgiin & Wasserman 2008 



Colaiuda et al. 2008) 



It is important to realise that upper limits on e cannot be used to constrain properties of QCD 
or the composition of the neutron star, which may simply have an ellipticity much lower than the 
theoretical maximum. The upper limits on e do, however, constrain the internal magnetic field to 
be less than of order 10^^ G, if Gas A is spinning fast enough to radiate gravitational waves in the 
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Rotation period / ms 



20.0 13.3 10.0 8.0 6.7 




Gravitational wave frequency / Hz 



Fig. 3. — Upper limits at 95% confidence (dots) on the intrinsic strain ho of gravitational waves 
from Cas A, and the indirect limit (line). The gravitational wave frequency is assumed to be twice 
the rotation frequency. Systematic uncertainties are not included; see Section [3] for discussion. 
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Rotation period / ms 
20.0 13.3 10.0 8.0 6.7 

5 .0 n ' ' ' ' I ' ' ' ' \ ' ' ' ' \ ' ' ' ' r 




100 150 200 250 300 

Gravitational wave frequency / Hz 



Fig. 4. — Upper limits at 95% confidence (dots) on the equatorial ellipticity e of Cas A, and 
the indirect limit (line). The gravitational wave frequency is assumed to be twice the rotation 
frequency. Systematic uncertainties are not included; see Section [3] for discussion. 
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Fig. 5. — Upper limits at 95% confidence (dots) on the amplitude a of r-mode oscillations of 
Cas A, and the indirect limit (line). The gravitational wave frequency is assumed to be 4/3 times 
the rotation frequency. Systematic uncertainties are not included; see Section [3] for discussion. 
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searched frequency band. 

The upper hmits on a (Figure [5| range from ~ 0.14 at 100 Hz (13 ms rotation period) 
to ~ 0.005 at 300 Hz (4.4 ms), assuming the canonical parameters of Eq. ([3]). If the r-mode 
amphtude varies with time, our hmits on a are rms values over the observing time and the indirect 
limits are rms values over the lifetime of the star. Our upper limits on a are within the range 
of runaway low-viscosity scenarios at all frequencies, and on the high end of the frequency band 
they are comparable to the finite- viscosity parametric instability thresholds which tend to serve as 
attractors for the evolution ( Bondarescu et ar]|2009 ). 

In several years the advanced LIGO and Virgo interferometers are expected to be in operation, 
with sensitivities an order of magnitude better than data searched here and extending to lower 
frequencies. Extrapolating from these results, a similar search on data from advanced interferom- 
eters would be expected to be sensitive to ellipticities of a few times 10~^, which are achievable by 
neutron stars without exotic matter or by stars with internal magnetic fields less than ~ 10^^ G, 
or r-mode amplitudes a few times 10~^. More sophisticated data analysis methods, such as hier- 
archical methods, would further increase the sensitivity. There are also more young non-pulsing 
neutron stars and other astrophysically interesting objects that could be targeted by a search of 
this type. 
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